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Abstract 



We obtain the renormalized equations of motion for matter and semi- 
classical gravity in an inhomogeneous space-time. We use the functional 
Schrodinger picture and a simple Gaussian approximation to analyze the time 
evolution of the X(p 4 model, and we establish the renormalizability of this non- 
perturbative approximation. We also show that the energy-momentum tensor 
in this approximation is finite once we consider the usual mass and coupling 
constant renormalizations, without the need of further geometrical counter- 
terms. 



Accepted to publish in Int. J. Mod. Phys.A 



1 



I. INTRODUCTION 



One of most rich and interesting subject in physics is the early universe, because it 
brings together two other topics that apparently have no connection: elementary particle 
physics and gravitation. According the standard cosmological model (Big Bang), the early 
universe was in such a high temperature state that particle physics must be used to describe 
the matter in this condition. The meaning of this is that the evolution of universe in the 
initial instants depended upon the particle spectrum and their interactions. In particular, 
the dynamics of scalar fields could give rise to the appearance of an inflationary era @ , that 
occurred during and/or after one of many phase transitions that the universe experienced. 

In general, the study of inflationary scenarios is made considering only some aspects of 
scalar field dynamics responsible for the existence of an inflationary phase. For instance, 
there are analyses of the classical effects of inhomogeneities in a scalar field in homoge- 
neous spaces 11, as well as in inhomogeneous spaces H. On the other hand, the studies of 
quantum and thermal effects do not that into account the consequences of departures from 
homogeneous situations @@. 

The purpose of this paper is to obtain the equations of motion for the A0 model and semi- 
classical gravity, taking into account non-perturbative quantum effects and inhomogeneities. 
This work is a generalization of ref. which analyzed the time evolution of homogeneous 
scalar fields in a Robertson- Walker space-time. This was accomplished in this reference using 
the functional Schrodinger picture and a simple Gaussian approximation. In this case it was 
possible to show the renormalizability of the equation of motion and energy-momentum 
tensor. In order to carry out a consistent analysis of the time evolution of inhomogeneities 
we must consider inhomogeneous spaces, however, we assumed that space is inhomogeneous 
in a single direction % to the problem be tractable. Furthermore, the exact solution of time 
evolution of an interacting system is impossible to be obtained, except for system whose 
Hamiltonian is quadratic, and consequently, approximation methods are needed. Therefore, 
we use variational approximations which lead to tractable equations to the various 
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parameters defining a Gaussian Ansatz. Although we work with a restrict class of states, 
the use of this variational methods allow us to obtain some non-perturbative effects through 
self-consistent equations. 

The inclusion of quantum effects in the evolution of the system gives rise to divergences 
in the equations of motion and the energy-momentum tensor. Therefore, we obtain the 
renormalization prescription for the mass and coupling constants in order to remove the 
infinities from the equations of motion. We also show that this renormalization prescrip- 
tion is enough to render the energy-momentum tensor finite in our approximation without 
the introduction of further geometrical counter-terms. In brief, our procedure leads to a 
consistent set of finite equations to matter and semi-classical gravitation in the presence of 
inhomogeneities. 

In this work we make use of the functional Schrodinger picture % to set up the equations 
of motion for this initial value problem. The conventional formulation of quantum field 
theory in terms of causal Green's functions in the Heisenberg picture is not specially suited 
to time-dependent problems that employ an initial condition for specific solution. Green's 
functions contain all information needed to determining transition rates, S"-matrix elements, 
etc., of system in equilibrium where initial data are superfluous. However, following a 
system's time evolution from a definite initial configuration is more efficiently accomplished 
in the Schrodinger picture description, where the initial data consist of specifying a state. 

For bosonic fields, the functional Schrodinger picture is a generalization from ordinary 
quantum mechanics to the infinite number of degrees of freedom that constitute a field. 
Therefore, we can use the mathematical/physical intuition acquired in quantum mechanics 
to analyze field-theoretic problems. However, the functional Schrodinger picture is not as 
widely used in actual calculations as the Green's function method since renormalization 
is more easily carried out in the latter framework. Nevertheless, it has been established 
renormalizability to the Schrodinger picture for both static and time- dependent cases. 
In this work we expand the list of situations where the renormalizability of this framework 
can be shown. 
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This work is organized as follows: In Section |I| we present the functional Schrodinger 
picture for a scalar field, as well as the variational principle used to obtain the equations 
of motion of the system. Section [TTT] contains the A0 4 model in which the homogeneous 
hypothesis in one direction (z direction) is abandoned and also our simplified Gaussian 
Ansatz. In section IV we renormalize the effective action in the Centrella- Wilson metric, 
which leads to the renormalized equations of motion for the variational parameters. The 
analysis of the renormalization of the energy-momentum tensor is shown in Sec. V. We draw 
our conclusions in Sec. VI. 



II. FUNCTIONAL SCHRODINGER PICTURE 

In the Schrodinger picture @, scalar field states are described by wave functionals \I/(0) 
of a c-number 0(x) at a fixed time, 

w) = (m ■ (i) 

The inner product is defined by functional integration, 

(#i|# 2 >4» jD<f>*l(<j>)* 2 (<t>), (2) 
while operators are represented by functional kernels 0(0,0'), 

0|*> <S> J D(j)' O(0,0')^(0') . (3) 

For the canonical field operator at a fixed time 0(x) (the time argument is common to all 
operators in the Scrodinger picture, so it is suppressed), we adopt a diagonal kernel $(x) =^> 
0(x)<5(0 — 0'); the canonical commutation relations determine the canonical momentum 
kernel to be II(x) =>- (l/i)[5/5<j)(x.)]5(<f) — 0'). In this way $ acts by multiplication on 
functionals of and II acts by functional differentiation. Hence, the action of any operator 
constructed from II and $ is 



1_5_ 

i 5d) 



9(4) ■ ( 4 ) 
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The fundamental dynamical equation is the time-dependent functional Schrodinger equa- 
tion for the time-dependent state functional ^(<f); t). This equation takes definite form, once 
a Hamiltonian operator H(U, $) is specified: 



In order to the initial value problem be completely defined, we must also specify the initial 
wave functional. 

In the case of a free scalar field it is rather simple make the connection with the usual 
approaches to field theory in curved space-time fl. For instance, the vacuum is given by a 
Gaussian wave functional which is not unique since there are many orthogonal solutions to 
the equations of motion. Moreover, we can also make field combination to define creation 
and annihilation operators to discuss particle production in curved space-time. For further 
details see refs. || , @ . 

The time-dependent Schrodinger, Eq. (|5|), cannot be solved exactly, unless the system 
is described by a quadratic Hamiltonian. To obtain a solution for a non-linear (interact- 
ing) system we shall employ a variational approximation in which Dirac's time- dependent 
variational principle @ is implemented approximately. 

In classical mechanics the Newtonian equations are derived by Hamiltonian's variational 
principle, which requires stationarizing the classical action I = J dt L(q,q), 



with the constraint that the variation 5q(t) vanishes at the end points of the integral defining 
the action. 

The quantum analogue of this principle for pure states is the Dirac's variational principle, 
which leads to the time-dependent Schrodinger equation. First, we define the effective action 
T as the time integral of the diagonal matrix element of id/dt — H: 




(5) 




(7) 
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and then we demand that T be stationary against arbitrary variations of and subject 
to appropriate boundary conditions, like it occurs in the Hamilton's variational principle. 
Applications of this variational principle with restricted variational Ansatz, in the Rayleigh- 
Ritz manner, result in tractable self-consistent dynamical equations for the parameters used 
in the Ansatz, which still retain some of the non-linearity of the complete problem. 



III. THE MODEL 



Many analyses of the time evolution of the scalar field associated to inflation were car- 
ried out assuming that the space is homogeneous and isotropic. In this work, we relax 
this hypothesis in order to assess the importance of inhomogeneities. To accomplish this 
goal and keep the problem tractable, we assumed that space is inhomogeneous only in one 
direction, which we defined to be the z direction. This resulting space-time is the sim- 
plest inhomogeneous one. This scenario is conveniently described by the Centrella- Wilson 
metric B, 

/ ?. , A?.nrr2 , a 9. n 7.2 a?, nr. n a ?. n7. \ 



V 



a 2 + A 2 (3 x2 + A 2 f3 z2 


A 2 f3 x 





A 2 (3 Z 


A 2 (3 X 


A 2 














AH 2 





A 2 (3 Z 








A 2 



(8) 



The parameters A, a, (3 X , f3 z , and I are functions of the time t and the coordinate z. The 
a term is a lapse function and it measures the orthogonal proper time interval aSt between 
neighboring slices a coordinate distance 5t apart. A shift vector, (3 l = ((3 x ,0,j3 z ), is used to 
simplify the three metric and keep it diagonal. The three metric of the slices is 7^ and is 
take to have the diagonal form 



^1 0^ 



in = A 2 



l 2 
1 



(9) 
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The function I is an anisotropic metric shear variable and the conformal factor A is a 
generalization of the usual scale factor R in the homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) models. Notice that the signature of the metric (§) is (-,+,+,+). 
The diagonal form of 7^ is preserved when the coordinate trajectories are adjusted by using 
/3 l , which results in two nonzero shift vector components in the four metric (||). 

The choice of this metric help us understand inhomogeneous, anisotropic cosmologies in 
terms of a simple, homogeneous, isotropic FRW model. For example, this metric generalizes 
the fiat (k = 0) FRW cosmologies, for which A(z,t) = R(t), a 2 = 1, (3 X = (3 Z = and 
1 = 1. As we can easily see, the number parameters defining this metric is substantially 
larger than the ones used to describe Robertson-Walker space-times, even assuming such a 
minimal modification. 

We assumed that the relevant degrees of freedom for the analyses of inflation are described 
by a A(/> 4 model, whose action is given by 



/ = J d n x C = - J d n x^~g r-g^d^dv® + + 



(10) 



whereas n = d + 1 with d being the number of spatial dimensions. From this action, we 
obtain that the canonical momentum conjugated to field $ is 

dC 



n = S = -v /= ^(^ $ + ^ $ • (H) 
Using the last two equations we can calculate the Hamiltonian density of our model: 

w(n,$) = 



°' IT 2 + <7y*n<9;$ + 7T ^- 7= (g ij + cr 2 g it g jt ) d&dfi 



2v^ 2a 
oJl V 2 + 4! 



(12) 



where e is AH 2 . 

We used the functional Schrodinger picture to describe the states and a time-dependent 
variational approximation based on the Gaussian Ansatz 



(0, t) = N(t) exp ji £ tt(x, t) [0(x) - y?(x, t)] 

xexp|-£ [0(x)-^(x,t)] [^- 1 (x,y,t)-^(x,y,t)] [0(y) - y?(y, t)]} , (13) 



where the variational parameters are <p, it, Q, and E and we abbreviated the integral in 
d spatial dimensions as J* x = / d d x. The physical meaning of the parameters of this wave 
functional can be inferred from linear and bilinear averages. The linear averages are given 
by 

($(x))=^(x,t) , (14) 

(II(x))=7r(x,t) , (15) 

while bilinear averages are 

($(x)$(y)> = y?(xMy) + fi(x, y, t) , (16) 
(n(x)n(y)) = 7r(x)7r(y) + V^x, y, t) + 4 (EQE) (x, y, t) , (17) 
($(x)n(y)) = ~5(x - y) + 2 (QE) (x, y, t) , (18) 

where we have used the matrix notation (OK,) (x, y) = J z C(x, z)/C(z, y). Moreover, from 
the average value of the operator i(d/dt), appearing in the effective action (0), we find that 
the imaginary part of the covariance function (E) plays the role of a canonical momentum 
conjugate to the real part Q, and the pair (tt, ip) has the same meaning, 

= jT7r(x,t)v9(x,t) + E(x,y,t)f2(y,x,t) . (19) 



Gaussian Ansatz in Curved Space 

In the spirit of Ref. ||, we simplified the Gaussian Ansatz specifying the momentum 
dependence of the kernels Q and E. This simplification is necessary not only to reduce the 
task of proving the renormalizability of our model, but also to make possible to apply it to 
numerical simulations. The kernels cam be expressed as a Fourier transform (TF) 

fi(x,y,t)= / e^<*-^tt(k,t) , (20) 

and we choose that 
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ft(k,t) 



2Je 



(21) 



E(k,t) 



c 



(22) 



8 {e [A- 2 (kl + /- 2 A;2 + A; 2 ) + a] } 

where the momentum-space integral (/ d d k/(2iT) d ) is denoted by / k . Therefore, we trade 
the variational parameters Q and E by a and (, in order to simplify our analysis. The 
parameter n is constrained by the renormalizability of the model, as we shall see below. It is 
interesting to notice that this form for the kernels corresponds to homogeneous fluctuations 
in the physical coordinates S. This is a reasonable hypothesis since the high frequency 
fluctuations should be insensitive to large scale effects. 

The Gaussian effective action is obtained through the substitution of the Ansatz given 
by Eqs. (|T3D, fl2~T|), and ( |22"D into the Dirac's variational principle (0) with the Hamiltonian 
density (0), 

1 - 2 



r = dt 

J Jx 

+ i(3 x ' 
+ if3 z 

+ h 



a 



7T(p 



2^1 

d^5 d {^-y) lx= y 



7r - /3 z nd z ip - oyfk ■ 



2A 2 



|x=y 



J {^(x, y, t) + 4*£(x, y, t)} c^(y, x, t) 
| {^(x, y, t) + 4*£(x, y, t)| ^(y, x, t) 



(£0) (x, x, t) - ^= | -IT 1 (x, x, t) + 2(SfiS) (x, x, t) | 

+ ^f{^ )2 + r ^ X)2 + ^ )2 }^ X ' y ' t ) 
fi(x, x, t)0(x, x, £) > . 



- ft 



2A 2 
2 AoVe 



x=y - {l? + ^ 2 J fl(x,X,t) 



(23) 



Notice that we have set ft = 1 throughout this paper except in (|23|), in order to understand 
the structure of the Gaussian approximation. Once more, the effective action have quantum 
corrections up to O (ft 2 ) , however, the variational kernels are non-trivial functions of ft, 
which re-introduces partially all orders in ft. In order to evaluate the Gaussian effective 
action we considered ip = tp(z,t), and used the following expectation values, 
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{(dj<$>) n) = -ndjtp + 2i {GdjQ) (x, x, t) , (24) 

(IldjS) = — ftdjtp - zof V(x - y) | x=y +2i (Gfyfi) (x, x, t) , (25) 

(fl^fl,*) = {djip) - ^ x) 9] y) fi(x, y, t) | x=y , (26) 
where we defined 

G(x,y,t) = ^- 1 (x,y,t) + 4zS(x,y,t) . (27) 



The stationary points of the effective action (f23|), with respect to the variational param- 
eters a, (, 7r, and (p, leads to the equations of the motion of these parameters. Notwith- 
standing, the Gaussian effective action exhibits infinities in the limit d — > 3 due to the short 
distance behavior of Q. Therefore, we must analyze its renormalization to obtain finite 
equations of motion. 

IV. RENORMALIZATION OF EFFECTIVE ACTION IN THE 
CENTRELLA- WILSON METRIC 



We regularized the divergences that appear in the effective action (|2"3j ) using dimen- 
sional regularization on the spatial dimension d, since this is the natural procedure in the 
Schrodinger picture @. In this section we show that this effective action become finite by 
the renormalization prescription H 

i = J_ I 1 

A X R (47r)( rf + 1 )/ 2 (d- l)(3-d) ' 1 J 

Initially we shall concentrate on the divergent terms of the effective action (^), which 

are the ones involving the kernels Q and E. We divided the divergent terms into two sets 

that are finite by themselves after the substitution of the renormalization prescription. The 

first set contains the divergent terms in the last three lines of the expression fl2~3|) , 

D\ = --Vr^x,*) + W 2 + rW 2 + d^ 2 } n(x, y ,t) 



>< 2 A 2 

A 

2* 



x=y 



f/i 2 + -if 2 ] Q(x, x, t) - ^^fi(x, x, t)ft(x, x, t) . (30) 
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In order to obtain the divergences in D±, we used the Fourier transform of the kernels, 
leading to 



a 



+ 



f I } \ 1 r (k 2 x + l~ 2 k 2 + k 



A- 2 (hi + l~ 2 k 2 y + k 2 ) + a 
1 



^ Jk J A- 2 (kl + l- 2 k 2 y + k 2 ) + a Jk ' J A- 2 (k' 2 + l~ 2 k' 2 + k' 2 



a 



> . 



(31) 



'A- 2 (kl + l~ 2 k 2 + k 2 ) + a 

The above expression is written in terms of co-moving coordinates. However, it is convenient 
to express it in the physical space variables, 



Px — k x /A , p y — ky/Al , p z — k z /A , 



(32) 



that leads to 



D 1 



cr\/e 



p 2 + a 



a 



1 2 X 2 



2 7 P Vp 2 + a 2 \ 2 I J P y/p r +~a 



X 

+ T77 



16 7 P ,p' y/p r +~a sjp 12 + a j 
Using dimensional regularization to evaluate the above expression we obtain 

A 



(33) 



awe 



a 



(4vr)( d + 1 )/ 2 (l -d) VA 2 
Xa 2 



2(4vr)( d + 1 )/ 2 (l -d) VA 2 / " V 2 

where A is an arbitrary mass scale which we choose to be given by 



(34) 



A 



2 _ ^ e -(7-V2) 

47T 



(35) 



Substituting the renormalization prescription ( p28|) - (29) into (J3J), we obtain Di in terms 
of renormalized quantities [x 2 R and Xr, which is finite in the limit d — > 3: 



1 I 2 64vr 2 



In 4 



2 , 2 



2X, 



2X 



R 



(36) 
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Notice that the above expression has a non-trivial dependence on \r due to the renormal- 
ization prescription that we used. 

The second set of potentially divergent terms appearing in ( p3|) is 

1 



D 2 = 



+ i(3 z 
2a 



di x) 5 d ( 



x- yj 



|x=y 



|x=y 



J {^(x, y, t) + 4*£(x, y, t)j d^Q(y, x, t) 
J y, t) + 4*£(x, y, t)} «9^(y, x, t) 



— (£0£)(x,x,t) + (SO) (x,x,t) . 
v e 



(37) 



The Fourier transform of fl37D , after the substitution of the Ansatze (|2l|) and ( 22|) can be 
written in terms of the physical momenta as 



c 



1 

i — — 



K 



2 16 e n {p 2 + ^ n+1/2 



a) 



((3 x Px + (3 z Pz ) 



32e r 

where we used that 
fi(x,x,t) = / Cl(k,t) 

Jk. 



2A 21 
a + a { — - + — 



64e 2n+3/2 (p2 + a ) 2 «+V2 
1 



A ' 3/ J J ( p 2 + ~,\ n + 3 / 2 



a 



fe 2A 2l\ 



(38) 



+ 



/2A 21 

a + a U + 37 



p (p 2 + a 



.3/2 



(39) 



In dimensional regularization integrals of the form J pi also vanish fl, while the well- 
known result 



-111 



nil p 



n 



(40) 



l P (p 2 + 2p-P~M 2 ) a (2n) d T(a)(-P 2 -M 2 ) a - n/2 V" 2, 

shows that the above integrals / Pij (p 2 + a) n+1/2 also vanish since they correspond to the 
case P = 0. Now the only potentially divergent terms in D 2 depend upon n. Since D 2 does 
not contain jj, or A, that could absorb a divergence, we choose n > 2 in £, in order to avoid 
new divergences in Z^- Consequently, (|38]) is finite in the limit d — » 3 and given by 



finite 



< 2 r ,_C_ 

64e 2 "+ 1 271+1/2 32e" 



/2A 2f 

d + a U + 37 



^n+3/2 + 



2 A 2f 
~A~3l 



'n+1/2 



(41) 
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where we defined 



'^iw^f- (42) 

Finally, substituting ( |36| ) and (^TJ) into ( |23| ) we obtain the renormalized effective action 
T R = jdt £ {vr0 - 01 + 0**0^ - ^ («9^) 2 + + D* | 



+ 



32e n 



d + a 



2* ° + 64^ 



'2A 21 
~A + 37 

1 



ta ^ 2 a ^ r 

^> ~ 64e 2n+l 7 2n+l/2 



/e 2i 2/'\ 



a 2 In 



1| («-/4) 2 ! /4 

2 J 2Ar 2A# 



(43) 



At this point we can obtain the finite (renormalized) equations of motion for the variational 
parameters n, (p, ( and a, by varying (^) with respect to them. 



5ip 
STr 



=> n = - 



ay/e axp , 



0^0 = (3 z d z ip - 



2^~e 



,-j I , 1 2A 2/\ _ 

=> a = J n+3/2 \ ^iW/J " " 1 T + 37 J ™ +3/2 



'e 2A 2l\ 
.e"T" 37 4+1/2 



(44) 
(45) 

(46) 



5Tr 

5 a 



C = CJ3/ 2 |C(" 4 

+ c(™ + 

- 32ae n+1 / 2 



d + a 
1 



e 

2 I \ e ~ 74 
2i 2Z 
T + 37 
1 



2 A 2f 



31, 

In+5/2 



n+3/2 



2^ 32vr 2 



a In [ 

V/4, 



A 



(47) 



V. RENORMALIZATION OF ENERGY-MOMENTUM TENSOR IN 

INHOMOGENEOUS SPACE 



The variation of the action fllPf ) with respect to the metric tensor g^, results in the 
energy-momentum tensor 
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(48) 



In this section we shall analyze the renormalization of the expectation value of T^ v in our 
Gaussian Ansatz. 

A. Renormalization of Diagonal Components of Energy-Momentum Tensor 



The expectation value of T M „, for fi = u, in the Gaussian state (|T3|) is given by 
- 2 2ap z 



7T 



(T a ) ='—+l3 z2 {d^f + ^7Td z¥ > 

e v e 



(-a 2 + A 2 j3 x2 + A 2 P 

P x \-idV>6 d {x - y) | x=y +42 (Gd x Q) (x, x, t) 



a 
a 



P z \-idP5 d (x - y) | x=y +4z (Gd z Q) (x, x, t) 



e )2 + r 2 5« 2 + e )2 l n(x,y,t) 



|x=y 



- ^-fi- 1 (x, x, t) - ^ (V + ^ 2 ^j fi(x, x, £) - ^(x, x, t)fi(x, x, £) | 

2 

+ ^( a 2 + A 2 l3 x2 + A 2 f3 z2 ){mZ){x,x,t) , 



(49) 



vr 2 /3* 



\2 , ^ 2 2 , ^ 



<?W = | ~ + ^ (d^) 2 + ^ 2 + y 



5 w 2 + r ^(*) 2 + 9 (x) 2 i n(Xj y; t) u=y + n -i( X) x , t) 



2A 2 
1 



+ | [ V 2 + ^(x, x, t) + ^n(x, x, t)fi(x, X, t) 



(SfiS)(x,x,t)l-9W 2 fi(x, y) t) 



x=y j 



(50) 



(Tyy) = -{Mf {-^ + |J (^ + ^ 2 + ^ 



2A 2 



9W 2 + l~ 2 d^ 2 + d^ 2 } Q(x, y, f) | x=y + g^T 1 (x, x, t) 
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+ \U 2 + ^vA ^(x,x,t) + ^n(x,x,t)fl(x,x,t) 



(£fi£)(x,x,t)} -r 2 ^ 2 fi(x,y,t) | x=y 



(51) 



(T«> = (^y?) 2 - A 2 
1 



7T 



(,9^) 2 + ^L 2 + * 4 



/3 2 



2A 2 



a w +r 2 9 (x) +9 w n(x ) y ) *)|^ r +-n- 1 (x,x,t) 



+ \U 2 + \v 2 J ^(x, x, t) + ^fi(x, x, *)fl(x, x, t) 



- ^(EOE)(x,x,t)}-aW 2 0(x,y,f) 



x=y 



(52) 



The divergent terms that are common to the four diagonal components are given by 



a w 2 + r a a« a + d ^) 2 ] n( x , y, t) | x=y --o-^x, x, t) 



~ \ V + 2"^ ) X ' ^ ~~ x ' *M X > x > *) 



(53) 



where there is no sum in the \i index. Substituting the Ansatz ([H|) and using the physical 
momenta results in 



D 



a 



1 ' 2 , \ 2 



A 



1 



1 



(54) 



32 7 P ' i/p' 2 + a j V P VP 2 + a 

At this point we evaluate this expression using dimensional regularization in order to extract 
the divergences, 

9 m (a\( d - 3)/2 -f3-d" 



D 



< !l < ~ (4tt)( c( + 1 )/ 2 (1 -d) \A 2 J 
Xa 2 



r 



.2 , \„2 



a ~ V + 2 V j a 



a Y d ~ 3)/2 T (3-d' 



(55) 



2(47r)( d + 1 )/ 2 (l -d) \A 2 J ~ V 2 

Notice that the ($ 4 ) term gives rise to a double pole since its expectation value is propor- 
tional to the square of fl(x, x, t). However, this double pole comes multiplied by A, which 
transform the double pole into a single one due to (|29|). 

Substituting the renormalization prescription ( |28|) - ( |29|) for /i 2 and A, we can see that 
the the only terms that survive in the limit d — > 3 are 
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D 



32^ 



16tt 2 



, of 



A 



(56) 



which still contains a pole. Nevertheless, there are other divergences in ( (4*9| ) - ( |52] ) that 
should be added to this one. In the case of (T u ), the remaining divergent terms are 

'2 



D« = - 



/Ftytf + p*ep*\ fi(x, y, o | x=y +^fi' 1 (x, x, t) 



(57) 



We already knows to calculate the terms in Q57p, adopting the same procedure used in the 
previous section, 



a 



-a z + A 2 (3 x2 + A 2 j3 z2 



In 



a 



2 I 2 



(58) 



32tt 2 V ^ ^ J [(3-d) 2 [ V^fi 

Summing (0), (|58j) , we can to see that the coefficient of the pole in the limit d — > 3 is 
canceled, when \i — v — t: 



D u + D 3 = (-a 2 + A 2 f3 x2 + A 2 (3 Z 
(a - fi 2 R f fj.% \ 



j a 2 


a 2 


In 




1" 




64tt 2 




~ 2 



+ 



(59) 



2A# 2X R J 

Here it is easy to verify that the remaining terms in fl4~9"|) are finite, which leads to the 
following for (T u ) in terms of renormalized coupling constants fi R and A#: 



(T tt ) R = - + (3 z2 (d zV f + ^*d a tp 



+ 



-a 2 + A 2 /T 2 + A 2 /T 2 

a 2 
64tt 2 

+ (a 2 + A 2 (3 x2 + A 2 /3 



m 4 Ui 



7T 2 1 o a > 

2^ 2 ,,4 ~1 



(a-/4) + /4 



2A 



R 



2A 



;2 



g4 e 2n+l 

The divergences in (T^-), when i = j, can be written as 



^2n+l 



/2 



(60) 



Ai - ^^ x)2 Q(x,y,i) | x=y , 



(61) 



where there is no sum in the % index. The Fourier transform of second term in ( |6~ID is given 
by 
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whose value in the limit d — > 3 is 

p2 



gu f V _ guor I i 1 
2d J P ~ 32vr 2 1 (3 - d) 2 



In 



a. 



2 I 2 



(63) 



In this way, the coefficient of the pole in (BTJ) vanishes, resulting in 



In 



D u - —d\ fl(x,y,t) | x=y = g u I --<p - — 



a 

,/4 



2 / 2 



+ 



{a - fi 2 R f fx 4 R 



(64) 



2Xr 2Xr t 

The spatial diagonal components of (T^) in terms of renormalized constants coupling fi 2 R 
and Xr, are given by 



+ 



a" 



64tt 2 



2e ' 2A 2V ~ Z ^ 64e 2n+i /2 ™+i/2+ 2 V 2 





i" 


Hi) 






~ 2_ 



(a -Mb) + Mr 



2A; 



2A 



R 



{Tyy)R — 1 2 {T XX ) 



R 



(T zz )r = {d z (p) 2 + {T xx )i 



(65) 
(66) 
(67) 



B. Renormalization of Off-Diagonal Components of Energy-Momentum Tensor 



First of all, we will analyze the off-diagonal elements of (T^), when one of their index is 
the time t. Considering the Gaussian state (|1|), these expectation values are 



7T 



(T xt ) = -'- + A 2 f3* 



1 fa \2 i M 2 2 i ^ 4 



2za 
7^ 



+ 



2A 2 
1 



(G^ft) (x, x, t) + A 2 (3 X ( - (EfiE) (x, x, t) - ift- 1 (x, x, t) 

I e 8 

+ r 2 5W 2 + d^ 2 - 2A 2 p*dWdW] n(x, y, f) | x=y 



> i M 2 + A(x,x,t) - ^fi(x,x,t)fi(x,x,t) 



x=y > 



(68) 
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(Tyt) 



2ia 



(G0 w fi)(x,x,t) - \p x d^d^+(3 z d^d^} fi(x,y,t) | x=y 



(69) 



(T zt ) = ^=nd zV -^ + A 2 p z 



2^ +Y V + 4! </9 



(G9 z n)(x,x,t) + a 2 /? 2 f-(£fi£)(x,x,t) - -n-^fox,*) 

I e 8 



+ 



2ia 
1 

2 ( ^ + 2^ ) fi ( x ' x ' *) ~ g fi ( x ' x > *M X > x , t) 



d^ 2 + r 2 <9« 2 + d^ 2 - 2A 2 /3^W^] fi(x, y, t) 



|x=y 



- ^w 2 n(x, y ,*) 

We use the expectation values 



|x=y 



(70) 



in order to obtain (T ty ) 



(71) 



(T ty ) = (T yt ) = ^=df ) 5 d (x-y) 

Following the procedure used in the previous subsections, we can write (|68| ) and (|70|) in 
terms of renormalized coupling constants jj? R and \r, 

r2 ( 



(^xt)ii = (T tx ) R = ~ + #F l~ (d zV f + e 



a 



6 4 e 2n+l /2 ™+ 1 /2 - 2^ Z 



a 



64tt 2 



{TztlR — (Ttz) 



a 



In 



Tr- 



et 

k/4 



2 i 2 



+ 



2A 



2A 



(72) 



a: 



R 



I <^ 2 

yjl * 2e \ 2A 2 v ^ ; 64e 2n+1 + / 2^ 

a - /4) 2 a4 1 



of 
64^ 





i" 








~ 2 



+ 



2A 



2A 



R 



(73) 



while (T y t) give us an integral that vanishes in the momentum space and, therefore, 

(T y t) R = . (74) 

The spatial components of off-diagonal elements of (T^J) are easily calculated: 

(Tij) = (Tji) = ^(x, y, t) | x=y . (75) 

It is easy to see that (T^) vanishes too since Q75|) is also a null integral in the momentum 
space, like as (ffll) . 



In the same way of ref . || , notice that the expectation of value (T^) become finite in our 
Gaussian approximation, once we take into account the mass and coupling renormalization, 
and that n > 2. This means that the infinities which exist in the free scalar field model are 
canceled by the presence of interactions, more specifically the term —g^ u (A/8) @. 

In addition, almost all divergences that appear when A# = (free theory) vanishes in 
consequence of the structure equation of motion for a. The only infinite that survives is 
guvlAtl^Ri which is absorbed by a redefinition of the cosmological constant. 



The functional Schrodinger picture is a useful tool to the analysis of initial value problems 
like the evolution of the early universe during and/or after an inflationary period. In this 
work we establish that a simple Gaussian variational approximation gives a consistent picture 
of the time evolution of matter and gravitation when the space is inhomogeneous in one 
dimension. In fact we were able to prove the renormalizability not only of the equations of 
motion of the scalar field, but also of the energy-momentum tensor. In order to obtain a 
finite effective action we used the same renormalization prescription that renders finite this 
quantity in the homogeneous case S. Moreover, the energy-momentum tensor is finite once 
we substitute the bare couplings and masses by the renormalized ones - that is, we do not 
need to introduce further geometrical counter-terms. This is a remarkable result taking into 
account that our approximation includes some non-linearities through the self-consistent 
equations of motion. 

In principle, the finite equations of motion for the scalar field supplemented by the semi- 
classical Einstein equation, 



VI. CONCLUSIONS AND DISCUSSION 



G 



^Gn{T^ v )gr , 



(76) 



where (T^ U ) GR is given by 



/4 



(77) 



2X R ' 
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may be used to study dynamical questions about the universe for inhomogeneous situations, 
such as stability of the space and the conditions for inflation setting in. Nevertheless, we 
should keep in mind that the Gaussian approximation in higher-dimensional field theory 
suffers from well known shortcomings, analogous to the that appear in the large- N approxi- 
mation. For instance, it does not describe spontaneous symmetry breaking correctly in four 
dimensions. In addition to that the behavior of the kernels assumed in the simplified Ansatz 
should describe well only the high frequency modes. Therefore, we are confined to the study 
of the chaotic inflation scenario within the framework developed here . 
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